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Abstract.  Starting  from  an  equation-of-motion  analysis  of  the  density  matrix,  phenomenological 
kinetic  equations  are  derived  that  allow  a study  of  semiclassical  intraband  and  quantum-mechanical 
tunneling  contributions  to  the  current  densities.  An  analytical  expression  is  obtained,  which  de- 
scribes the  field  dependence  of  the  electron  densities  in  the  upper  and  lower  minibands.  It  is  shown 
that  a population  inversion  does  not  occur  as  long  as  the  scattering  time  for  thermal  generation 
is  larger  than  the  interminiband  relaxation  time.  Numerical  results  for  the  current  density  are 
discussed. 


1.  Introduction 

Recently,  there  have  been  intensive  investigations  of  the  influence  of  external  electric  fields 
on  semiconductor  superlattices  (SLs).  The  physics  of  biased  SLs  is  extremely  rich  due  to 
the  large  number  of  parameters  that  can  be  controlled  quite  independently.  Most  of  the 
interest  has  been  caused  by  the  need  of  understanding  the  interplay  of  Bloch  oscillations, 
Zener  tunneling  between  Bloch  bands,  and  interaction  effects.  From  a theoretical  point 
of  view  the  simplest  models  involve  only  a single  miniband.  This  approach  leaves  out  all 
interesting  interminiband  phenomena,  which  are  captured  most  simply  in  a two-band  tight- 
binding  model.  The  purpose  of  the  present  work  is  to  study  the  nonequilibrium  carrier 
kinetics  by  an  equation-of-motion  analysis  of  the  density  matrix  and  a derivation  of  an 
expression  for  the  current,  from  which  semiclassical  intraband  and  quantum-mechanical 
tunneling  contributions  can  be  identified.  Starting  from  a microscopic  theory,  we  derive 
phenomenological  kinetic  equations  by  using  the  constant  relaxation-time  approximation. 
The  dynamics  of  electrons  is  treated  in  a simple  two-band  SL  model.  Former  theoretical 
studies  [ . ] of  unipolar  devices  will  be  extended  by  taking  into  account  all  relevant 
horizontal  and  vertical  electron  transitions.  In  addition,  the  widely  used  Esaki-Tsu  model 
[ 1 is  extended  to  two-band  systems. 

2.  Theoretical  model 

The  quantum-kinetic  equation  for  the  Wigner  transformed  elements  of  the  density  matrix 
/,!’  (k)  (with  v,  i>'  being  miniband  indices),  whose  explicit  spatial  dependence  describing 
field  domains  is  not  taken  into  account,  is  given  by  [ J 

j}Eyk  + -£,.(fe)]J  fv'(k)  (1) 

E [e^(fe)/;'(fc)  - e».v(fc)/,r(fc)]  = E E 

v /'A  fci 
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Fig.  2.  The  dynamical  Kerr  rotation  spectra  as  a function  of  time  delay  in  the  femtosecond  time 
domain. 


Fig.  3.  The  spectral  dependence  of  the  induced  Kerr  rotation  for  two  delay  times  t = 0 and 
1 = 500  fs. 
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where  the  matrix  elements 


0w,'(fc)  = + K)Wk:x*(k  + K)  (2) 

K 

detemiine  the  wavefunction  overlap  calculated  from  the  SL  envelope  functions  Xu-  The 
0-term  in  Eq.  ( 1 ) describes  electric-field-induced  vertical  transitions  via  Zener  tunneling. 
We  consider  a two-band  model  characterized  by  the  tight-binding  minibands 

trk2,  A i 

£i(k)  - — ^ -I-  — (1  -cosk:d).  (3) 

2 m 2 

t)2k2  At 

si(k)  — \ 1 + fp  + A | + — ^(1  +cos  k-d).  (4) 

2m  2 

Ai  ( A 2 ) is  the  bandwidth  of  the  lower  (upper)  miniband,  ands^  the  energy  gap  between  the 
minibands.  k±  is  the  quasi-momentum  of  the  lateral  electron  motion.  To  lowest  order  in  the 
coupling  constant,  only  the  transition  probabilities  W,1,1.  w£,  Wjl,  Wj,1.  W\\.  and  W jy 
have  to  be  accounted  for.  We  will  study  the  kinetic  equation  ( 1 ) in  the  constant  relaxation 
time  approximation  by  making  use  of  the  replacement 

ReW|?(fc',  k)  = ReWf;(k’.  k ) -»  -8k'k/ r (5) 

for  the  interminiband  transition  probabilities.  The  horizontal  interminiband  transitions 
via  thermic  generation  and  recombination  W£,1  are  described  by  relaxation  times 
t[  and  respectively.  Finally,  for  the  description  of  intraminiband  transitions,  we  use 
the  extended  Esaki-Tsu  model  f ].  Focusing  on  high  electric  fields,  where  £2  > l/ren 
(£2  = eEd/h  is  the  Bloch  frequency  and  reff  an  effective  scattering  time),  and  taking  into 
account  the  periodicity  of  /,!'  (k)  along  the  field  direction  k- . we  switch  to  the  Wannier- 
Stark  representation  by  Fourier  transforming  the  distribution  function  | ] 

GO 

//(fc)  = e'^/.r'ffex./).  (6) 

/=  — OO 


We  neglect  the  smooth  k-  dependence  in  Ql_ll,’(kz  ) and  obtain  the  following  analytical  result 
for  the  electron  density  «2  of  the  upper  miniband 

m - + T/Ti  A _ y-'  <?(02  7 

"2  4(£2t)2A0  + r/r,' + r/rr  ° ^ (/£2r  - QgT)2  + 1 


In  the  expression  for  Aq,  we  introduced  the  abbreviations 


<?(/)  = (— i i QnJi 


/A,  + A2\ 
V 2/i£2  ) 


- 


Eg  + (A]  -t-  A2)/2 


— £2(022  — 0n)-  (8) 


with  J/  denoting  the  Bessel  function.  In  the  absence  of  thermal  generation  r[  —*■  oo  and 
under  the  condition  /t2  <§;  1.  we  recover  the  result  n2  % 2£22tttAq  derived  already  in  [ ]. 
From  Eq.  (7),  it  is  seen  that  population  inversion  occurs  only  under  the  unrealistic  condition 
Tt  > t{  independently  of  the  value  of  the  electric  field  strength. 
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Fig.  1.  (a)  Electric-field  dependence  of  the  carrier  density  for  the  upper  (solid  lines)  and  lower 
(dashed  lines)  minibands.  The  bandwidth  of  the  lower  (upper)  miniband  is  A|  =5  meV  (A:  = 
40  meV ).  The  two  minibands  are  separated  by  an  energy  gap  of  = 100  meV  and  the  temperature 
is  T = 300  K.  We  use  identical  intraminiband  relaxation  times  ri  = r?  = 0. 1 ps  and  the  following 
values  for  the  interminiband  scattering  times  r — 0.3  ps  and  r{  = 20  ps.  For  the  solid  curves  from 
bottom  to  top.  the  interminiband  recombination  time  r,  is  given  by  1 . 3,  and  5 ps.  The  positions  of 
Zener  resonances  at  IQ  — Qv  (I  — 1. 2)  are  marked  by  thin  vertical  lines. 

(b)  The  same  as  in  (a)  for  the  relative  current  density  ./- /./o  with  /u  = ens(  A|  + At  ) /2 h.  ns  is  the 
carrier  sheet  density. 


An  analytical  expression  for  the  current  density  is  obtained  from  [ ) 

jz  = Ij  £i  <*:)/'  (k’)WH  (V.  k)  + (Mk:)  fhk’)WH(k\  k) 

kk'  1 

+f  i (k:)  fhk')W;\  (ft\  k)  + e2(k-J  /,'  (fc')  W'fffc'.  fc)  J.  (9) 

where  n denotes  the  total  electron  density. 

3.  Numerical  results  and  discussion 

Numerical  results  for  the  carrier  densities  and  the  relative  current  jz/h  ate  shown  in  Fig.  1 (a) 
and  (b).  The  set  of  parameters  used  in  the  calculation  is  given  in  the  caption.  The  matrix 
elements  Q^,,'  have  been  estimated  with  wavefunctions  of  a quantum  well  with  infinite 
potential  barriers.  A characteristic  feature  in  these  graphs  is  the  appearance  of  Zener  reso- 
nances at  IQ  = Qg.  which  are  shifted  to  larger  field  strengths  with  increasing  energy  gape^ 
so  that  at  low  electric  fields  the  well-known  Esaki-Tsu  current-voltage  characteristics  [ 1 
is  recovered.  The  lineshape  of  the  Zener  resonances  depends  sensitively  on  the  miniband 
widths.  Whereas  the  carrier  density  n 2 is  a function  of  the  sum  A|  + A2,  pronounced 
current  resonances  occur  only  for  narrow  lower  and  broad  upper  minibands.  With  increas- 
ing interminiband  relaxation  time  r the  peaks  become  sharper.  When  the  generation  and 
recombination  lifetimes  are  identical  (rj  = r-J).  we  obtain  from  Eq.  (7)  n 2 — 0.5  inde- 
pendently of  the  electric  field  strength.  As  long  as  thermal  generation  can  be  neglected 
(ti  < Tj ),  a population  inversion  cannot  occur  in  the  considered  SL  model. 
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